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Abstract. Recent cosmic microwave background (CMB) 
measurements over a large range of angular scales have be- 
come sensitive enough to provide interesting constraints 
on cosmological parameters within a restricted class of 
models. We use the CMB measurements to study inflation- 
based, cold dark matter (CDM) critical density universes. 
We explore the 4-dimensional parameter space having 
as free parameters, Hubble's constant Ho, baryonic frac- 
tion fift, the spectral slope of scalar perturbations n 
and the power spectrum quadrupole normalization Q. 
We calculate minimization values and likelihood in- 
tervals for these parameters. Within the models consid- 
ered, a low value for the Hubble constant is preferred: 
Ho = 30^7* km s~"'^Mpc~"'^. The baryonic fraction is not 
as well-constrained by the CMB data: = O.OTtom- The 
power spectrum slope is n = 0.91^q^2- The power spec- 
trum normalization is Q = 18 ± 2.5 /iK. The error bars 
on each parameter are approximately la and are for the 
case where the other 3 parameters have been marginal- 
ized. If we condition on n = 1 we obtain the normalization 
Q = 17 ± 1.0 /j,K. The permitted regions of the 4-D pa- 
rameter space are presented in a series of 2-D projections. 
In the context of the CDM critical density universes con- 
sidered here, current CMB data favor a low value for the 
Hubble constant. Such low- Ho models are consistent with 
Big Bang nucleosynthesis, cluster baryonic fractions, the 
large-scale distribution of galaxies and the ages of globular 
clusters; although in disagreement with direct determina- 
tions of the Hubble constant. 

Key words: cosmic microwave background — cosmology: 
observations; theory 
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1. Introduction 

The standard picture of structure formation relies on the 
gravitational amplification of initially small perturbations 
in the matter distribution. The origin of these fluctua- 
tions is unclear, but a popular assumption is that these 
fluctuations originate in the very early universe during an 
inflationary epoch. The most straightforward incarnation 
of this inflationary scenario predicts that the fluctuations 
are adiabatic, Gaussian, Harrison- Zeldovich (n = 1) and 
that the Universe is spatially flat (Kolb & Turner 1990). 
To avoid violating primordial nucleosynthesis constraints, 
the Universe should be dominated by non-baryonic mat- 
ter. The cold dark matter (CDM) model has been the 
preferred model in the inflationary scenario (Peebles 1982, 
Liddle and Lyth, 1993). 

The statistical properties of CMB fluctuations provide 
an ideal tool for testing CDM models. CMB data offer 
valuable information not only on the scenario of the ori- 
gin of cosmic structures, but also on the early physics of 
the Universe and the cosmological parameters that char- 
acterize the Universe. Using the CMB to determine these 
parameters is the beginning of a new era in cosmology. 
This truly cosmological method probes scales much larger 
and epochs much earlier (z > 1000) than more traditional 
techniques which rely on supernovae, galaxies, galaxy clus- 
ters and other low-redshift objects. The CMB probes the 
entire observable universe. 

Acoustic oscillations of the baryon-photon fluid at 
recombination produce peaks and valleys in the CMB 
power spectrum at sub-degree angular scales. Measure- 
ments of these mo del- dependent peaks and valleys have 
the potential to determine many important cosmolog- 
ical parameters to the few percent level (Jungman et 
al. 1996, Zaldarriaga et al. 1997). Within the next decade, 
increasingly accurate sub-degree scale CMB observations 
from the ground, from balloons and particularly from two 
new satehites ( MAP: Wright et al. 1996, Planck Sur- 
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veyor: Bersanelli et al. 1996) will tell us the ultimate 

fate of the Universe (fio), what the Universe is made 
of {fib, ^cdm) and the age and size of the Universe 
{h = -ffo/(100/ km s~^Mpc~^) with unprecedented pre- 
cision. 

In preparation for the increasingly fruitful harvests of 
data, it is important to determine what the combined 
CMB data can already tell us about the cosmological pa- 
rameters. In Lineweaver et al. (1997), (henceforth "paper 
I"), we compared the most recent CMB data to predic- 
tions of COBE normalized flat universes with Harrison- 
Zel'dovich (n = 1) power spectra. We used predominantly 
goodness-of-flt statistics to locate the regions of param- 
eter space preferred by the CMB data. We explored the 
h — fib plane and the ft, — Aq plane. 

In the present paper we focus on the range for h fa- 
vored by the CMB in the context of CDM critical density 
universes and wc broaden the scope of our exploration to 
the 4-dimensional parameter space: h, fib, n and Q. Our 
motivation for choosing this 4-D subspace of the higher 
dimensional parameter space is that (i) it is the largest 
dimensional subspace that we can explore at a reasonable 
resolution with the means available and (ii) it is centered 
on the simplest CDM model: fig = I, = 0, n = 1. This 
model is arguably the simplest scenario for the formation 
of large-scale structure. One of our goals is to see what 
is required of such a model if it is to explain the current 
set of large-scale structure data, and what could eventu- 
ally force us to accept the fine-tuning demanded by the 
inclusion of another cosmological parameter, such as the 
cosmological constant. 

Bubble's constant Ho is possibly the most important 
parameter in cosmology, giving the expansion rate, age 
and size of the Universe. Recent, direct, low-redshift mea- 
surements fall in the range [45-90] but may be subject 
to unidentified systematic errors. Thus it is important 
to have different methods which may not be subject to 
the same systematics. For example, CMB determinations 
of h are distance-ladder-independent. Current CMB data 
are not of high enough quality to draw definitive model- 
independent conclusions, however in the restricted class of 
models considered here, the CMB data are already able to 
provide interesting constraints. 

The quantity fib is important because we would like 
to know what the universe is made of and how much nor- 
mal baryonic matter exists in it. The combination flbh'^ is 
relatively well-constrained by the observations and the- 
ory of primordial nucleosynthesis, but the uncertainty 
on the Hubble constant means that the value of fib is 
rather poorly constrained. The question of just how many 
baryons there are in the Universe has received close at- 
tention recently due to estimates of the baryon fraction 
in galaxy clusters and attempts to constrain fib by mea- 
suring the deuterium in high-redshift quasar absorption 
systems. 



The parameter n is the primordial power spectrum 

slope that remains equal to its primordial value at the 
largest scales (low i). It is important because it's mea- 
surement is a glimpse at the primordial universe. Although 
generic inflation predicts n = 1, a larger set of plausible in- 
flationary models is consistent with 0.7 ^ n ^ 1.0. Model 
power spectra and particularly the amplitude of the first 
peak depend strongly on n. Thus, an important limita- 
tion of paper I was the restriction to n = 1. By adding 
n as a free parameter we obtain observational limits on n 
and quantify the reduced constraining ability of the CMB 
observations when n is marginalized. 

The power spectrum quadrupole normalization Q is 
important because it normalizes all models. Here we treat 
Q as a free parameter. 

We examine how the contraints on any one of these 
parameters change as we condition on and marginalize 
over the other parameters. We obtain minimization 
values and likelihood intervals for h, fib, n and Q. As in 
paper I, wc take advantage of the recently available fast 
Boltzmann code (Seljak and Zaldarriga 1996) to make a 
detailed exploration of parameter space. 

All the results reported here were obtained under a 
restrictive set of assumptions. We assume inflation-based 
CDM models of structure formation with Gaussian adia- 
batic initial conditions in critical density universes {fig = 
1) with no cosmological constant (Ao = 0). We have ig- 
nored the possibility of early reionization and any gravity 
wave contribution to the spectra. Wc do not test topo- 
logical defect models. We use no hot dark matter. We 
have used the helium fraction Yhc = 0.24 and a mean 
CMB tcmpcrturc To = 2.728 K. Although wc have not 
yet looked carefully at how dependent our results are on 
these assumptions we make some informed estimates in 
Sect. 6.1 where wc also discuss previous work using similar 
data sets and similar methods to look at different families 
of models. 

In Sect. 2 wc describe the data analysis. In Sect. 3 we 
present our results for h and fib and discuss their depen- 
dence on some plausible variations in the data analysis. 
We discuss non-CMB constraints and compare them with 
our results in Sect. 4. In Sect. 5 we present our results for 
n and Q. In Sect. 6 we add some caveats and summarize. 

2. Method 

2.1. Data 

The data used are described in paper I, however we have 
updated some data points and now include several more 

measurements : 

• updated Tenerife point (Hancock et al. 1997): ST^ff = 

34.1+1^^3^ at 4// = 20, 

• added BAM point (Tucker et al. 1997): STeff = 
55.6tlli at 4// = 74, 

• updated the two Python points (Piatt et al. 1997): 
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STeff = GOZll at ^eff = 87 and STeff = GGZH at 
4// = 170, 

• added MSAM single- and double-difference points 
(Cheng et al. 1996) STeff = 40.71?^;^ at 4// = 159 

Snjj = UAtlli at 4// = 263. 

With the exception of the Saskatoon points, the cali- 
bration uncertainties of the experiments were added in 
quadrature to the error bars on the points. The MSAM 
values are the weighted averages of the first and second 
MSAM flights. The error bars assigned encompass the la- 
limits from both flights. In paper I. wc did not include the 
MSAM points because of possible correlations with the 
Saskatoon results. However the MSAM results are sub- 
stantially lower than the Saskatoon results in this crucial 
high-£ region of the power spectrum and it is not clear that 
avoiding such correlations is more important than the ad- 
ditional information provided by the MSAM data. The 
figures presented in this paper include the MSAM points 
however we have also performed calculations without 
MSAM. We discuss these results in Sect. 3.1. 

2.2. Calculation 

A two-dimensional version of our calculation is de- 
scribed in paper I. In this work we generalize to 5 di- 
mensions and use a likelihood approach to determine the 
parameter ranges. We treat the correlated calibration un- 
certainty of the 5 Saskatoon points as a nuisance param- 
eter "u". For each point in 5-D space we obtain a value 
for x^(/i, r^b, n, Q, zi). We assume u comes from a Gaus- 
sion distribution about its nominal value with a disper- 
sion of 14%. This Gaussian assumption amounts to adding 
[(1 — u)/0.14]^ to the calculation described in paper 
I. For example, in paper I, our notation Sk-14, SkO and 
Sk-|-14 corresponds to u = (0.86,1.00,1.14) respectively 
(however we did not assume a Gaussian distribution and 
so did not add the extra factor to the x^ values). 

For each point in the 4-D space of interesting param- 
eters, u takes on the value which minimizes the (Avni 
1976, Wright 1994). At the minimum in 4-D, Xmi„, the pa- 
rameter values are the best-fit parameters. To obtain error 
bars on these values, we determine the 4-D surfaces which 
satisfy x'^ih,^b,n,Q) = Xmin + ^' where x = [1, 4, 9, 16]. 
Under the assumption that the errors on the data points 
are Gaussian (cf. de Bernardis et al. 1996), the x = 1 el- 
lipsoid can be projected onto any of the axes to get the Icr 
confidence interval for the parameter of that axis. If one is 
not interested in 1-D intervals but rather in the confidence 
regions for 2 parameters simultaneously, then one would 
use X = [2.3, 6.17, 11.8, 19.3] (see Press et al. 1989 for 
details). To make the figures, we project the 4-D surfaces 
onto the two dimensions of our choice and obtain contours 
which we project onto either axis. 

We have normalized the power spectra using the con- 
version Cii = Q^^^, where is the power spectrum 
output of the Boltzmann code. 
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Fig. 1. Likelihood contours in the h — Clb plane from re- 
cent CMB measurements. We have conditioned on n = 1 and 

Q = 17 ^K. The 'X' marks the minimum. The contours are 
at levels Xmin + ^ where x = [1, 4, 9, 16] (see Section 2.2). 
When projected onto either of the axes these regions give the 
approximate size of the 1, 2, 3 and 4a confidence intervals re- 
spectively. The area within the la contour has been shaded 
dark grey. The resuh, h = 0.351!^;^^, = 0.02^°;°*, is given 
in Table I. Big Bang nucleosynthesis predictions favor the light 
grey band defined by 0.010 < Qbh^ < 0.026 (see Sect. 4.1). 



3. Results for h and fit 

The permitted regions of the 4-D parameter space are pre- 
sented in a series of 2-D projections which contain like- 
lihood contours from a combination of the most recent 
CMB measurements. There are four groups of figures cor- 
responding to the four planes h — Qb, h — n, h — Q and 
n — Q] Figures 1 - 3, 5, 6 - 7 and 8-9 respectively. The 
thick 'X' in each figure marks the minimum. Areas within 
the Ict contours have been shaded. 

The best-fit values and confidence intervals displayed in 

the figures are summarized in Table I which thus con- 
tains the main results of this work. The values of h, fi;,, 
n and Q at the minimum of the 4-D x? are given with 
error bars from the projection onto 1-D of the Xmin + t 
surface. In Figures 1 through 4 the region preferred by 
Big Bang nucleosynthesis (BBN) is shaded light grey 
(0.010 < n^h? < 0.026, see Sect. 4.1). 
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.20 .30 .40 .50 .60 .70 .80 .90 1.0 

h 

Fig. 2. Contours and notation as in previous figure except 
here the normahzation Q has become a free parameter, i.e., 
for a given h and f^b, Q takes on the vahic that minimizes 
the value of at that point. The minimum and the la errors 
on h and f2(, are the same as in Figure 1 however the 2, 3, 
and 4a contours are noticeably larger. The higher values of Q.h 
permitted here correspond to Q < 17 /uK. 

In Figure 1 we have conditioned on n = 1 and Q = 
17 AtK. The results arc h = 0.35to°l, = 0.02lg;g^. At 
2(7, 0.25 <h< 0.60 and 0.0 < ^t, < 0.13. 

The contours and notation of Figure 2 are the same 
as in Figure 1 except that we have let the normalization 
<5 be a free parameter. That is, for each value of h and 
Oft, Q takes on the value (within the discretely sampled 
range) which minimizes x^{h-, ^b)- The minimum and the 
1(7 errors on h and Ob are the same as in Figure 1; h stays 
low. The 2, 3, and Aa contours are noticeably larger than 
in Figure 1. Within the 2a contours, the higher values of 
Oft permitted correspond to Q < 17 /xK. 

Figure 3 displays the main result for h of this paper. 
The result is more general than the results of Figures 1 
and 2 since no restrictions on n and Q are used. Exam- 
ining Figures 1, 2 and 3 sequentially, the dark grey la 
contours can be seen to get larger as we first condition on 
and then marginalize over n and Q. With both n and Q 
marginalized we obtain h = 0.30tg J^ and Oft = 0.07^^:o7 
where the error bars are approximately la. At the min- 
imum, n = 0.91 and Q = 18 jjK. The value at this 
minimum is 16. The number of degrees of freedom is 23 
(= 27 data points - 1 nuisance parameter - 3 marginalized 
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.20 .30 .40 .50 .60 .70 .80 .90 1.0 

h 

Fig. 3. Same as previous figure except here both Q and 
n are free parameters. The preferred value of h stays low: 
h = O.SOIq J^ but at 2a, h can assume all values tested, i.e., 
the projection of the 2a region onto the h axis covers the entire 
axis. We also find Q.b = 0.071° q*. At the minimum n = 0.91 
and Q = 18 /iK. In contrast to the previous figure the higher 
values of Q.b permitted correspond to Q > 17 jxK and n < 1. On 
the right, the high h, Q.^ ~ 0.15 region within the 2a contour 
has n ~ 0.7 and Q ~ 20 /xK. 



parameters). The probability of finding a value this low 
or lower is ^ 15%. Thus the fit obtained is "good". 

The region of the /i — Oft plane acceptable to both 
the BBN and CMB have low values of h. Large values 
of h, especially in the BBN region (lower right of plot) 
are not favored by the CMB data. However at 2a, h is 
unconstrained since the projection of the 2a contours onto 
the h axis spans the entire axis. 

Since n is a free parameter the amplitude, but not 
the location, of the Doppler peak can vary substantially. 
This explains the shape of the 68% region allowed by the 
data: the possible range of O5 is enlarged, but not the 
confidence interval of the Hubble constant, which seems 
to be predominantly determined by the position of the 
Doppler peak. 

The disjoint 2a contour region on the right is char- 
acterized by the parameter values h <i 0.65, Oft ~ 0.15, 
n f» 0.7 ± 0.1 and Q w 20 ± 2 /jK. The minimum of this 
region is at /i = 1.0. 
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Fig. 4. This plot has no CMB information in it. The bands 
are the constraints from 4 non-CMB measurements discussed 

in Sect, 4: Big Bang nucleosynthesis ("BBN"), cluster bary- 
onic fraction ("Clusters"), the galaxy and cluster scale density 
fluctuation shape parameter ("F") and the age of the oldest 
stars in globular clusters ("Age": region between the verti- 
cal lines). To quantify the combination of these 4 constraints, 
we show contours obtained by assuming for each constraint 
a two-tailed Gaussian probability around the central values. 
The 1 and 2a contours are thus the result of an approxi- 
mate joint likelihood of 4 non-CMB constraints for Qo ~ 1, 
Ao = CDM models, A low h is preferred. This was pointed 
out in Bartlett et al. (1995), The combination of these four, 
independent, non-CMB measurements yields ~ 40 ± 7 and 
Qb « 0,18 ±0,06, This region of the h — Qh plane is very similar 
to the region preferred by the CMB data (compare Figure 3), 
The point of this diagram is not to show that h is low since we 
have ignored the numerous, more direct, z ^ measurements 
of h which find h w 0,70 ± 0,10, The point of this diagram is 
to show that an important set of independent constraints also 
favors low values of 7i in fio = 1, Ao = CDM models. 



out the MSAM data points. When the MSAM points are 

not included the h and fi;, minimum and contours do not 
change significantly. For example the results from Figure 
3 without MSAM are h = O.SOtom and fife = O.OStg og. 
• Saskatoon calibration treatment 

We have treated the Saskatoon calibration as a nuisance 
parameter from a Gaussian distribution around the nom- 
inal Saskatoon calibration with a dispersion of 14%, The 
values of u at the minimum in this technique are 
u ^ 0.82. Netterfield et al. (1997) have compared the 
Saskatoon results to the MSAM first fiight results in the 
north polar region observed by both experiments. They 
find a best-fit calibration of —18% which is equivalent to 
the 0.82 discussed above. Thus, there is some evidence for 
a lower nominal Saskatoon calibration 

However, preliminary results based on a relative cali- 
bration between Jupiter and Cas A at 32 GHz imply that 
a +5% ± 7% Saskatoon calibration is appropriate (Leitch 
et al. 1997). Using this calibration changes the results 
slightly. For example, the analog of Figure 3 yields tighter 
contours around the unchanged h minimum: h = O.SO^Q oy 
and Qfe = 0.19 with la error bars larger than the range 
probed. The preference for fi;, ^ 0.15 — 0.20 is increased 
(independent of h) and the avoidance of the high h, low 
Ofe, BBN region is increased. The value of the minimum 
increases from « 16 to « 21 thus the fit is still good; 42% 
probability of having a lower . 

We have also let m be a free parameter from a uni- 
form distribution, i.e., a free-floating Saskatoon calibra- 
tion. The analog of Figure 3 yields h = 0.30^0 07' = 
O.Obtom- At the minimum u = 0.78(= -22%), = 14 
and the probability of obtaining a lower x^ is 7%. We 
have also obtained results assuming three different Saska- 
toon calibrations; the nominal value, 14% higher and 14% 
lower. The minimum x^ stays at /i « 0.30 in all three 
cases. 

Thus several plausible choices of data selection and 

data analysis producing ~ 20% variations in the ampli- 
tude of the Doppler peak, do not strongly affect the low 
h results from the GMB. The many measurements on the 
Doppler slope in the interval [10 < £ < 200] contribute 
strongly to determining the position of the peak and thus 
to the preference for low h (see Figures 3 and 4 of paper 
1). h Ki 0.30 appears to be a fairly robust CMB result for 
the critical density CDM models tested here. 



3.1. Robustness of results to data analysis choices 

Since there is a ~ 2(t inconsistency between the MSAM 
and the Saskatoon data and since there may be unknown 
systematic errors, we have performed some checks to see 
how dependent our results are on the various ways of an- 
alyzing the data. 
• Without MSAM 

The figures and Table I results include the MSAM data 
points, but we have also performed x^ calculations with- 



3.2. Ofe results 

Our CMB constraints on flf, are weaker than our con- 
straints on h; the la contours in Fig 3 are elongated verti- 
cally and yield fib = 0.07t.o_or- Comparing Figures 1 and 
2 with 3, one sees that it is the marginalization over n 
which opens the ili, ^ 0.07 region (where n < 0.9). White 
et al. (1996) highlight the merits of high baryonic fraction 
(fife ~ 10% - 15%) models. We confirm that the CMB 
X^ Ic region is centered near this range but the valley of 
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Table I: Parameter Results 
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" parameter values at the minimum. The error bars come 
from the projected Icr contours in the figures. 

" "free" means that the parameters were free to take on any 
value within the discretely sampled range which minimized the 
value of x^- As a result of the discretization, these reported 
minima can be displaced from the true minima by up to half a 
grid point on both axes, i.e.. Ho = 30 should be taken to mean, 
27.5 < Ho < 32.5. The underlying matrix of model points is 
described by 0.15 <h< 1.00, step size: 0.05, 0.01 < fib < 0.4, 
step size: 0.012, 0.5 < n < 1.5, step size: 0.06, 10 < Q < 30, 
step size: 1. Thus we have tested over 200,000 (18 x 34 x 18 x 20) 
models. 

Result from joint likelihood of non-CMB constraints (see Fig- 
ure 4). 



minima is very shallow. In the context of our models, non- 
CMB data can still constrain slightly better than the 
CMB. See our discussion of Figure 4 in Sect. 4.5 where we 
report ^6 » 0.18 ± 0.06. 



most published results. These limits are plotted in Fig- 
ures 1-4 and are labelled "BEN" in Figure 4. We use 
f2(, = 0.015 as a central value. The BBN constraints are 
independent of Oq; Aq, n and Q and thus do not depend 
on our = 1, Ao = assumptions. Lyman limit sys- 
tems yield a somewhat model-dependent lower limit for 
the baryonic density, lending support to the higher values 
of Ofe/i^ (Weinberg et al. 1997, Bi & Davidsen 1997). 

4-2. X-ray cluster baryonic mass fraction 

Observations of the X-ray luminosity and the angular size 
of galaxy clusters can be combined to constrain the quan- 
tity n;^^^^- We adopt the range 0.04 < ^^h^^^ < 0.10 
(White et al. 1993) with a central value of 0.06 (Evrard 
1997). These limits seem to be inconsistent with BBN if 
f^o = 1 and h ^ 0.50. This is known as the baryon catas- 
trophy and has led some to believe that < 1- The 
severity of this catastrophy in Qq = 1 models can be ex- 
amined in Figure 4 by comparing the "BBN" region with 
the "Clusters" region. For Clb ^ 0.15, Oo = 1 models allow 
consistency between the nucleosynthetic and cluster data 
for low values of h. 



4.3. Matter power spectrum shape parameter T 

Peacock & Dodds (1994) made an empirical fit to the mat- 
ter power spectrum using a shape parameter P. For fig < 1 
models, P can be written as (Sugiyama 1995) 



T = hflo exp 



Oo + l 



(1) 



We adopt the 2a limits of the empirical fit of Peacock & 
Dodds (1994)(see also Liddle et al. 1996a) and include the 
n dependence 



0.222 < P 



0.32 ( - 



1 < 0.293. 



(2) 



Under the assumption that 0.8 < n < 1.2, Equation 2 
becomes 0.169 < F < 0.373. This is the F constraint used 
in Figure 4. We use F = 0.25 as a central value. 



4. Non-CMB constraints in the h — fib plane 

Four independent non-CMB cosmological measurements 
constrain the acceptable regions of the h — ilb plane. 

4.1. Nucleosynthesis 

Primordial nucleosynthesis gives us limits on the bary- 
onic density of the Universe. Although deuterium mea- 
surements seem to be currently the most accurate bary- 
ometer, there is an active debate about whether they 
yield high (Tytler et al. 1996, Tytler & Buries 1997) or 
low (Songaila et al. 1994, Carswell et al. 1994, Rugers & 
Hogan 1996) baryonic densities. We have adopted the 
range 0.010 < Vti, < 0.026 because it encompasses 



4.4. Limits on the age of the Universe from the oldest stars 
in globular clusters 

There is general agreement that the Universe should be 
older than the oldest stars in our Galaxy. Thus a lower 
bound on the age of the Universe comes from an age de- 
termination of the oldest stars in the most metal-poor (= 
oldest) globular clusters. A reasonably representative sam- 
ple of globular cluster ages, tgc, in the literature is. 



Bolte & Hogan (1995) 


13.7 


< tgc 


< 


17.9 


Sarajedini et aJ. (1995) 


13.7 


< tgc 


< 


14.3 


Chaboyer (1995) 


11 


< tgc 


< 


21 


Jimenez ct al. (1996) 


11.5 


< tgc 


< 


15.5 


Salaris et al. (1997) 


10.4 


< tgc 


< 


14.0 
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Allowing ~ 1 Gyr for the formation of globular clusters, 

we adopt the range 11 < io < 18 Gyr with a central value 
of 14 Gyr. We use this relatively large interval to avoid 
overconstraining the models and to encompass most pub- 
lished results. Age determinations are Qq and Aq indepen- 
dent but converting them to limits on Hubble's param- 
eter depends on our fio = 1 and Aq = assumptions. 
In the models we arc considering here, our age limits are 
converted directly into limits on Hubble's constant using 
h = 6.52 Gyr/to which yields 0.59 > h> 0.36 (with the 
central value 14 Gyr corresponding to h = 0.47). This 
region is marked "Age" in Figure 4. 

4-5. Summary of non-CMB constraints and comparison to 

CMB constraints 

The constraints we adopt from BBN, cluster baryonic frac- 
tion, r, and the ages of the oldest stars in globular clusters 
(as described above) are. 



BBN 


0.010 


< 


Of, 


< 


0.026 


Clusters 


0.04 


< 




< 


0.10 


r 


0.169 


< 


r 


< 


0.373 


GC Ages 


0.36 


< 


h 


< 


0.59 



Bands illustrating these constraints are plotted in Figure 
4. Since these constraints are independent, it is not obvi- 
ous that they should be consistent with each other. They 
are consistent in the sense that there is a region of overlap. 
This consistency is improved if ^lbh^ turns out to be high 
as indicated by Tytler & Buries (1997). To visualize more 
quantitatively the combination of these four constraints, 
for each constraint we assume a two-tailed Gaussian distri- 
bution around the central values. This allows the flexibility 
to account for asymmetric error bars. We then calculate a 
joint likelihood, 



C{h, Qb) = CbBN * Cclus * Ct * Cage 

where, 

(Qb h? - 0.015)2 



CBBN{hMb) oc exp 
Cciua{h,^b) oc exp 
jCr{h,Clb) oc exp 
Cage{h) oc exp 



2 rr^ 

^ "BBN 



iSlb - 0.06)2 



2cr2 

clus 



{r{h,nb) - 0.25)2 



2gI 



{h - 0.47)2 



2<e 



(3) 

(4) 

(5) 

(6) 

(7) 
(8) 



fraction, F and stellar ages yield h w 0.40 ± 0.07 and 
rih « 0.18 ± 0.06 for fio = 1, Ao = 0. 

The point of Figure 4 is not to show that h is low since 
we have of course ignored the numerous, more direct, 2: ~ 
measurements of h which find h « 0.70 ± 0.10 (see for 
example Freedman 1997). The point of this diagram is to 
show that an important set of independent constraints do 
overlap and are consistent with each other for low values 
of h in the = 1, Aq = models considered here. 

Figure 4 should be compared with Figure 3 which has 

. There is an 



2 



The upper and lower limits of the four constraints de- 
termine the cr's for the two-tailed Gaussians. For ex- 
ample (TBBN,up = (0.026 - 0.015) and (JBBN,down = 

(0.015 - 0.010). The joint likchhood of the four terms 

C-BBN * Ccius * Ct * Cage is shown in Figure 4. The con- from Figure 3 since we are projecting the same 4-D 



contours of CcMB{h,Q.b) oc exp 
interesting consistency between the non-CMB constraints 
and the CMB constraints of Figure 3. Although they ex- 
tend over relatively small regions in the h — Qb plane, the 
1(7 regions of the non-CMB joint likelihood and the CMB 
overlap. 

In Figure 4 we see that the combination of four inde- 
pendent cosmological measurements indicate that a low 
value of h could make the CDM theory viable, as Bartlett 
et al. (1995) argued. The amplitude of small scale mat- 
ter fluctuations is an additional consistency check on this 
model. The value of ag at the minimum in Figure 3 is 
as = 0.54. This agrees quite well with values inferred 
from X-ray cluster data (Viana & Liddle 1996, Oukbir 
et al. 1996). 

Liddle ct al. (1996b) studied critical density CDM 
models. Based on the COBE normalization, peculiar ve- 
locity flows, the galaxy correlation function, abundances 
of galaxy clusters, quasars and damped Lyman alpha sys- 
tems, they found that h < 0.50 and n < 1 is preferred. 
Adams et al. (1996) come to similar conclusions. 

Bartlett et al. (1995) listed the advantages of a low h in 
critical density universes, the main point being that there 
exists a region of parameter space in which this simplest of 
models remains consistent with observations of the large- 
scale structure of the Universe. For example, there is the 
question of the age crisis. Recent h measurements point to 
values in the interval 0.60 < h < 0.80. In a critical density 
universe h = 0.70 implies an age of 9.3 Gyr, younger than 
the estimated age of many globular clusters, h = 0.30 
yields an age for the universe of 21.7 Gyr comfortably in 
accord with globular cluster ages. 

What is perhaps surprising is the fact that the CMB 
data do not rule out such a such a low value of Ho but 
seem to favor it within the context of this type of scenario. 
Of course, these low Ho values are in disagreement with 
current measurements of the Hubble constant. 

4.6. Other projections 

Figures 5 through 9 show our 4-D ellipsoid projected 
on to 2-D planes orthogonal to the h — flb plane. The limits 

obtained on h from Figure 5 are the same as we obtained 



tour levels are C{h,Qb)/Cmax = ea;j)(— 5[1,4]). The com- 
bined non-CMB constraints from BBN, cluster baryonic 



surface, h and n are positively correlated for h ^ 0.50 and 
possibly negatively correlated for h ^ 0.50. In Figure 6 we 
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Fig. 5. 

Q free. 

has Q ' 



Likelihood contours in the h~ n plane for fib free and 
h = O.SOl'o and n = O.Qlt",'^"- The h ^ 0.80 region 
- 20 ^K. At 2a, 0.55 < n < 1.23, and the largest n is 
for 0.50 < h < 70. In the la contour there is a strong h and n 
correlation: the lower values of h go with the lower values of n. 
When different values of n are allowed, the amplitude of the 
Doppler peak varies up and down but the location of the peak 
does not. 



see that with n fixed at 1, a high precision determination 
of Q is possible, Q = 17 ± 1.0 /uK. In Figure 7 the n = 1 
constraint is dropped. 



5. Results for n and Q 

Figures 8 and 9 show the strong anti-correlation between n 
and Q. This has been observed and discussed by many au- 
thors (e.g. Smoot et al. 1992, Seljak & Bertschinger 1993, 
Lineweaver 1994). This anti-correlation has a simple ex- 
planation. In the Ci{£) plot, any increase of the slope low- 
ers the y-intercept (at £ = 2) and any decrease of the 
slope raises the y-intercept. The anti-correlation is thus 
inherent with the use of the amplitude at = 2 as the 
normalization. 

The increase of the size of the error bars on n and Q as 
h and are conditioned on and then marginalized can be 
seen by comparing Figures 8 and 9. Our n and Q results 
are n = OM^oil and Q = 18±2.5AtK from Figure 9 where 
both h and have been marginalized. Conditioning on 
n = 1, we get Q = 17 ± 1.0 /uK (see also Fig. 6). 



Fig. 6. Likelihood contours in the h — Q plane for n = 1 and 
^6 free, h = 0.351!;;"?. Q = 17 ± 1.0 /xK. Notice that since n 
has been fixed at 1, the value of Q docs not rise in the high h 
region as it does in the next figure. The la contours are round: 
there is no correlation between h and Q. 



The four year COBE-DMR constraints on the ampli- 
tude and slope of the power spectrum at large angular 
scales are udmr = 1-2 ± 0.3 and Q = 15.3^2^ A*K, and 
conditioning on udmr = 1, Q = 18 ± 1.6 iJ.K (Bennett et 
al. 1996). These DMR results are in the context of = 1, 
Ao = CDM models and they need to be corrected due to 
the mildly model-dependent extended tails of the Dopplor 
peak which reach even into the low-^ region. After the 
correction, the DMR result becomes n « 1.05 ± 0.3 and 
Q « U.btli A«K. Thus our results from a combination of 
recent CMB measurements in the context of critical den- 
sity universes agrees well with the DMR-only result and 
reduces the error bars on both n and Q. 

For the standard CDM model of Figure 8 {h = 0.50, 

flo = I and fib « 0.05), we obtain 0.90 <n< 1.17 (95% 
CL) and 13 < Q < 20.5 /xK (95% CL). Using simi- 
lar methods and a similar data set, several authors have 
reported similar results, de Bernardis ct al. (1996) find 
1.0 < n < 1.26 (95% CL). White et al. (1996) find 
0.86 < n < 1.10 (95% CL). Hancock et al. (1997) find 
1.0 < n < 1.2 (68% CL). The variations of these stan- 
dard model n determinations are probably due to slightly 
different data sets, different treatments of the Saskatoon 
calibration and the details of the calculation. 
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.20 ,30 .40 .50 .60 .70 .80 .90 1.0 

h 

Fig. 7. Likelihood contours in the h — Q plane for n free and 
fif, free. Notice that in the high h region, Q ~ 20 /iK is pre- 
ferred. A comparison with Figure 5 shows that n and Q arc 
anti-correlated but this is seen more easily in the next two 
figures. 

6. Caveats cind summEiry 

6.1. Caveats 

Our low h measurements are in disagreement with current 
2 ~ measurements of the Hubble constant. Possible ex- 
planations for this discrepancy could be unidentified sys- 
tematics in the CMB data or the local h measurements. 
Galactic foregrounds could be a problem for the CMB 
while the distance ladder may need some readjusting for 
the local h measurements (e.g., Feast & Catchpole 1997). 
The best way to address these problems is with more and 
better data. This is being done rapidly. New detectors and 
better designed observations are improving the quality of 
both the CMB data and the more direct h measurements. 
SZ and supernovae observations are also increasing the 
redshift of the h measurements. 

A simple answer to the discrepancy between direct h 
measurements and our results is that the Universe is not 
well-described by the models considered here. It is possi- 
ble that one or more of our basic assumptions is wrong, 
or we could simply be looking at too restricted a region of 
parameter space. The shape of the primordial power spec- 
trum may be more complicated than the family of mod- 
els we are using. Inflation may be wrong and structure 



301 r 




n 

Fig. 8. Likelihood contours in the n — Q plane for h = 0.50 
and Qb = 0.05. Result n = 1.03l";!!g and Q = lltl'^o, or at 
the 95% CL, 0.90 < n < 1.17 and 13 < Q < 20.5 fiK. This 
standard model has been looked at by several other workers 
who obtain similar results (see Sect. 5). 



may not have formed from adiabatic curvature fluctua- 
tions. Topological defects may be the origin of structure. 

Wc emphasize that we have only considcrcid a partic- 
ular type of cosmological scenario, although arguably the 
simplest; the results we have presented here are valid un- 
der the assumption of inflation-based, Gaussian adiabatic 
initial conditions in a critical density universe (fio = 1) 
with no cosmological constant. We have not considered 
any early reionization scenarios or gravitational wave con- 
tribution. We have also not included any hot dark matter. 

Our fig = 1 assumption can be considered very restric- 
tive since plausible values for Oq in the range [0.2, 1.0] can 

change the power spectrum signiflcantly. For example, the 
position of the Doppler peak, ipeak is roughly proportional 

— 1/2 

to Qo ■ Thus low flo, by pushing Ipeak higher may per- 
mit higher h values to accomodate the high ipeak ^ 270 
of the CMB data. We are in the process of checking this 
assumption; we consider open models in Lineweaver et 
al. (1997, in preparation) 

Reionization models can affect the power spectra sig- 
nificantly by lowering the Doppler peak but this can be 
compensated by n values larger than 1. For example, de- 
Bernardis et al. (1996) have looked at reionization models 
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n 

Fig. 9. Likelihood contours in the n — Q plane for h free and 
Qfe free, n = OMt^ f^, Q = 18 ± 2.5 pK. The minimum is at 
h = 0.30, fib = 0.07, and u = 0.82. 



and find a best-fit early reionization at Zreion ~ 20 and 
n w 1.2. 

In paper I we took a brief look at flat-A models. The 
supernovae results of Perlmutter et al. (1997) can con- 
strain A better than the CMB data. The combination of 
the CMB, BEN and supernovae constraints in flat-A mod- 
els yields 0.23 <h< 0.72. 

If gravitational waves or any other effect plays an im- 
portant role in CMB anisotropy formation, wc expect that 
the inclusion of this effect in the family of models tested, 
will improve the resulting fits. However the inclusion of 
gravitational waves seems to make the fits slightly worse 
without changing the location of the minimum (Liddle 
et al. 1996b). Bond and Jaffe (1997) analyzed the com- 
bined DMR (Bennett ct al. 1996), South Pole (Gimdcrson 
et al. 1995) and Saskatoon (Netterfield et al. 1997) data 
using signal-to- noise eigenmodes. They looked at the pa- 
rameters h, n and erg in a variety of models. The inclusion 
of tensor modes for n < 1 models seems to produce small 
shifts in the likelihood surfaces. 

There may be extra-relativistic degrees of freedom (hot 
dark matter (HDM) or mixed dark matter (MDM)). dc- 
Bernardis et al. (1996) found that current CMB anisotro- 
phy measurements cannot distinguish between CDM and 
MDM models. We agree with this assessment and add 



that HDM and CDM models also cannot usefully be dis- 
tinguished with current CMB data. 

In addition to the h, fib, n and Q considered here, 
regions of a larger dimensional parameter space deserve 
further investigation including f2o, Aq, Qhdm ■ early reion- 
ization parameters such as zreion > tensor mode parameters 
ut and T, iso-curvature or adiabatic initial conditions and 
topological defect models with their additional parameters 
such as the coherence length. The fact that we obtain ac- 
ceptable values in our small 4-D parameter space lends 
some support to the idea that we may be close to the 
right model. However establishing error bars on broad- 
band power estimates is a relatively new science. If the 
Universe is not well-described by these models then as the 
data improve, work like this will show poor fits and 
other regions of parameter space may be preferred. 

6.2. Summary 

CMB measurements have become sensitive enough to con- 
strain cosmological parameters in a restricted class of 
models. The results we have presented here are valid under 
the assumption of Gaussian adiabatic initial conditions in 
a critical density (f2o = 1) universe with no cosmological 
constant. We have explored the 4-dimensional parameter 
space of h, Qb, n and Q. Our CMB-derived constraints on 
h, fib, n and Q exclude large regions of parameter space. 
We obtain a low value for Bubble's constant: Ho = SOl^^. 
The CMB data constrain only weakly: = Q.Qlto^^. 
For the slope and normalization of the power spectrum 
we obtain n = 0.91^°;^° and Q = 18 ± 2.5 /iK. The error 
bars on each parameter are for the case where the other 
3 parameters have been marginalized. When we condition 
on n = 1 we obtain the normalization Q = 17 ± 1.0 /xK. 

CMB constraints are independent of other cosmologi- 
cal tests of these parameters and are thus particularly im- 
portant. The fact that reasonable values are obtained 
means that the current CMB data are consistent with 
inflationary-based CDM models with a low Hubble con- 
stant. In the context of the models considered, the CMB 
results are consistent with four other independent cosmo- 
logical measurements but are in disagreement with more 
direct measurements of h. 
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